Spin dynamics in a strongly driven system: very slow Rabi oscillations 
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We consider joint effects of tunneling and spin-orbit coupling on driven by electric field spin 
dynamics in a double quantum dot with a multi-level resonance scenario. We demonstrate that 
tunneling plays the crucial role in the formation of the Rabi-like spin-flip transitions. In contrast to 
the linear behavior for weak electric fields, the spin fiip rate becomes much smaller than expected 
for the two-level model and shows oscillating dependence on the driving field amplitude in stronger 
fields. In addition, the full spin flip is very difficult to achieve in a multi-level resonant system. These 
two effects have a similarity with the Zeno effect of slowing down the dynamics of an observable by 
its measurement. As a result, spin manipulation by electric field becomes much less efficient than 
expected. 
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I. INTRODUCTION 

Fast reliable spin manipulation in quantum dots is one 
of the challenges in spintronics and semiconductor-based 
quantum information. The design of corresponding gates 
can be based on electric dipole spin resonance where the 
spin-orbit coupling (SOC) [I-lj] allows on-chip spin ma- 
nipulation by external electric field as well as electric 
read-out of spin states. [5| Without external driving, SOC 
effects on localized in quantum dots electrons are very 
weak and lead to long spin relaxation times, [y, 01 The 
important questions here are how fast can the gate oper- 
ate, what limits the manipulation rate, and how efficient 
is the spin manipulation in terms of the achievable spin 
configurations. [8| It seems that a stronger driving field al- 
lows for a faster spin manipulation, as predicted in a sim- 
ple Rabi picture of the driven oscillations. This picture is 
applicable for a single quantum dot with a parabolic con- 
finement, where the electron displacement from the equi- 
librium is linear in the applied electric field. [9| However, 
double quantum dots where tunneling plays the crucial 
role for the orbital dynamics, and the corresponding en- 
ergy scales are different from a single quantum dot, are 
more promising for observation of new physics and ap- 
plications in quantum information technologies, [l^ The 
tunneling makes the description of the SOC puzzling 
since the electron momentum is not a well-defined quan- 
tity at under-the-barrier motion, and the tunneling rate 
can become strongly spin-dependent. pll Il2| In addition, 
the double dots provide a possibility to study free and 
driven interacting qubits.flj, [Ij] Here we concentrate on 
one-dimensional systems attractive for spintronics [la - 
IVq and building quantum dots |18| - |20| and consider spin 
manipulation in single-electron double quantum dot |2l| - 
[26j by periodic electric field. 

We show that even for a basic quantum system such 
a single electron spin, the efficiency and time scale of 
the manipulation strongly depend on the electron orbital 
motion and, as result, to an unexpected dependence on 



the external electric field. [27| The nonlinearity of the spin 
and charge dynamics is expected to lead to unusual con- 
sequences on the driven spin behavior. [28| In a multilevel 
system Rabi spin oscillations are slowed down if the field 
is sufficiently strong, which challenges efficient spin ma- 
nipulation. We restrict ourselves to the single electron 
dynamics to demonstrate in the most direct way the non- 
trivial mutual effect of coordinate and spin motion on the 
Rabi oscillations. The slowing of the oscillations down 
at high electric fields is a truly unexpected general fea- 
ture of a multi-level system compared to the conventional 
two-level model and thus can occur in a broad variety of 
structures. 

This paper is organized as follows. In Section II we 
introduce quantum mechanical description of electron in 
a double quantum dot with spin-orbit coupling and mag- 
netic field. Section HI presents the model of driven dy- 
namics. In Section IV we apply the stroboscopic Flo- 
quet approach for the long-time evolution and obtain the 
properties of Rabi oscillations under various conditions. 
Conclusions of this work are given in Section V. 



II. MODEL, HAMILTONIAN, AND 
OBSERVABLES 

The unperturbed Hamiltonian Hq = k"^ /2m -f U{x) 
describes electron in a double quantum dot with the po- 
tential (see FigH]) Hi]: 



U{x) = Uo 



-'1 



(1) 



where k = ~id/dx is the momentum operator and h = 
I. The minima at — d and d arc separated by a barrier 
of the height Uq. In the absence of external fields and 
SOC the ground state is split into the doublet of even 
(tpg) and odd {^u) states. The tunneling energy AEg <^ 
Uq determines the tunneling time Ttun = 2tt/ AEg. The 
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FIG. 1: (Color online) (a) Qualitative picture of the spin 
dynamics induced by the interminima tunneling. (b,c) Free 
evolution of coordinate (solid line) and spin components (p-^) 
(dashed line), (uj,) (dashed-dot line) for B = 1.73 T (b) and 
B = 6.93 T (c). The initial state is the even combination of 
the states corresponding to the tunneling-split doublet with 
SOC taken into account. 



Zeeman coupling to magnetic field Hz 
\A.z\ is the Zeeman splitting. 
The SOC has the form: 



^zCTzf^, where 



ffs, 



{aDCTx + aRGy) k, 



(2) 



where the bulk-originated Dresselhaus (au) and 
structure-related Rashba (oifi) parameters determine the 
strength of SOC. In the presence of SOC the spatial par- 
ity of the eigenstates is approximate rather than exact 
being the qualitative feature of the coupling linear in the 
odd fc-operator, eventually resulting in the ability of spin 
manipulation by electric field. 

The quantities we are interested in are the spin com- 
ponents: 



{'^^it)) 



(Ji[x,t)dx, 



(3) 



where ai{x, t) = '4>^{x, t)aiXp{x, t) is the spin density, and 
expectation value of the coordinate 



(x(i)) = / xp{x,t)dx, 



(4) 



where p{x,t) = ip''{x,t)if}{x,t), and xp{x^t) is the two- 
component electron wavefunction. 

For numerical studies we diagonalize exactly the time- 
independent Hamiltonian Hq + Hso+Hz in the truncated 



spinor basis 'ipnix)\a), where ipn{x) are the eigenfunctions 
of iJo in Eq.([T]), and \a) where cr = ±1 corresponds to the 
spin parallel (antiparallel) to the z-axis, find correspond- 
ing eigenvalues, and obtain the new basis set IV'n)- We 
consider below as an example a GaAs-based structure, 
where the effective mass is 0.067 of the free electron mass, 
with d = 25-\/2 nm and Uq = 10 meV. In the absence of 
magnetic field the ground state energy is Ei = 3.938 
meV, and the tunneling splitting AEg ~ 0.092 meV, cor- 
responding to the transition frequency close to 23 GHz. 
To illustrate the spin dynamics, we consider a moderate 
external magnetic field with Az ~ AEg/2 correspond- 
ing to B = 1.73 T, and a relatively strong magnetic field 
with Az = 2AEg {B = 6.93 T) with the Lande factor 
g = —0.45. The parameters of the SO coupling are as- 
sumed to be an = 1.0 • 10~^ eVcm, and an — 0.3 • 10~^ 
eVcm, however, our results can be applied to various dou- 
ble quantum dots with different SOC parameters and 
thus have a quite general character. In particular, the 
change in the interdot barrier shape and geometry would 
modify only quantitatively the system parameters, in- 
cluding the energy levels, spinor wavefunctions, and, as 
a result, the resonant driving frequency. The increase in 
the interdot distance would decrease the tunneling split- 
ting, making such a system more sensitive to external 
influence from phonons, fiuctuations in the driving field, 
etc. 



III. DRIVEN DYNAMICS 

To demonstrate a nontrivial interplay of the tun- 
neling and spin dynamics, we begin with the coordi- 
nate and spin evolution of the electron initially local- 
ized near the —d minimum. Spin evolution of the state 
(iV'g) + IV'u)) |1) /V2 can be described approximately an- 
alytically taking into account four spin-split lowest levels 
and a simpler SOC Hamiltonian aRUyk as 



{<7x{t)) = OLR 



KAEg 



■ s\i\{A+t) s\\i{A-t) 



(5) 



where K = ~i{ipu \k\ ipg), which in the AEg ^ Uq limit 
can be accurately approximated as K = mdAEg, and 

El/ 4 + a%K^, where E± = AEg ± Az. Nu- 



A. 



merical results for coordinate and spin are shown in FiglT] 
With the increase in magnetic field, the effect of SOC 
decreases, leading to smaller amplitudes of precession, as 
can be seen from comparison of Figlljb) and Figdjc). In 
addition, both the initial state and spin precession axis 
change leading to a different phase shift between the ob- 
served spin components. 

Next we consider a periodic perturbation by electric 
field at i > 0: 



£{t) = £asm{u!zt). 



(6) 



Here ujz is the exact, taking into account SOC, frequency 
of the spin-fiip transition. For the chosen set of param- 



eters u)z is very close to A^. The field strength is char- 
acterized by parameter / defined as e£o = / x Uo/2d, 
where e is the fundamental charge. For the chosen sys- 
tem parameters, / = 1 corresponds to the electric field of 
approximately 1.5 x 10"^ V/cm, similar to Ref.Q- Here 
we consider different regimes of the strength and see how 
the change in the shape of the quartic potential produced 
by the field becomes crucially important for the spin dy- 
namics in two sets of energy levels produced by magnetic 
field. We build in the obtained |-0„) basis the matrix of 
the Hamiltonian V = cxEq smiujzt) and study the full 
dynamics with the wavef unctions: 



t/;(x,i) = ^^nWe 



-?:B„t 



IV-n) 



The time dependence of (,nit) is given by: 



dt 



ze£{t)J2^iit) 



•^In^ 



-i(Ei^E^)t 



(7) 



(8) 



where xin = {'4'i\x\ipn)- There are two different types 
of xin'- (1) matrix elements of the order of d due to the 
different parity of the wavefunctions in the absence of 
SOC, and (2) those determined by the SOC strength. 
In the weak SOC limit |Az - A^g| > aR\K\, the SO- 
determined matrix element of coordinate in the lowest 
spin-split doublet can be evaluated as 



2dKa 



h? 
'-^z 



i^E.y 



(9) 



In our calculations we assume that the initial state is 
the ground state of the full Hamiltonian, that is ^i (0) = 1 
and ^n>i(0) = 0. The entire driven motion of the system 
can be approximately characterized as a superposition of 
two types of transitions: resonant "spin-flip" transitions 
with the matrix element of coordinate determined by the 
SOC and off-resonant "spin-conserving" transitions with 
a larger matrix element of coordinate. J30j Both types are 
crucially important for the understanding of the spin dy- 
namics. With the estimate K sa mAEgd, in both cases 
considered by us {Az = Eg/2 and Az ~ 2Eg), we ob- 
tain d « lOxso- As a result, the off-resonant transitions 
are not weak compared to the required once. Through- 
out the calculation we neglect orbital and spin relaxation 
processes assuming that the driving force is sufficiently 
strong to prevent the decoherence on the time scale of 
the spin spin-flip transition. It is known that the peri- 
odic field forms a well-established driven dynamics even 
in the presence of damping as long as the level structure 
is not deeply disturbed by the broadening. For our pa- 
rameters it means that one can expect the observation 
of the predicted results in the currently available semi- 
conductor structures at temperatures moderately below 

iK.fg 

We begin with presentation of the short-time dynamics 
of coordinate {x)/d and spin (ax) for four initial periods 
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FIG. 2: (Color online) (a) Scheme of the levels and spin com- 
ponents involved in the Zeeman resonance, B = 1.73 T (left) 
and B = 6.93 T (right), (b) short-term time-dependent {x)/d 
and {(Jx), as marked near the plots, induced by external field 
with / = 0.1, B=1.73 T, Tz=88 ps; (c) same as in (b) for 
B=6.93 T, rz=22 ps. 



of the driving field (Figl2|). These resuts were obtained 
by the explicit numerical integration of Eq.(I5]) with a 
time step on the order of lO^^Tz- The other component 
(az) changes much slower and will be treated later on a 
long timescalc, which is the primary topic of our inter- 
est. It can be seen in Figl2]that the fast oscillations are 
accompanying mainly the local variations of observables, 
especially of the spin. Considerable changes such as 
Rabi oscillations of spin can be achieved only after many 
periods of the driving field. We will focus on this slow 
dynamics below. 



IV. FLOQUET STROBOSCOPIC APPROACH 

To consider the long-term time dependence of the pe- 
riodically driven system we apply the Floquet approach 
B [31n23 in the stroboscopic form. Here we remind 
the reader main features of this approach developed in 
Ref . [3a| . As the first step, the one-period propagator 
matrix \Jin{Tz) is obtained by a high-precision numerical 
integration of the system ^ at one period of the driv- 
ing Tz = 2tt/ljz in the basis of all unperturbed states. 



For numerically accurate Ui„(Tz), we obtain its eigen- 
values Eq which are the quasienergies of the driven sys- 
tem, and the corresponding orthogonal eigenvectors Af. 
As a result, the one-period propagator XJiniTz) can be 
presented as: 



V,,iTz) = Y.Af{AQye^ 



iEnTz 



(10) 



Its A'^-th power obtained by taking into account the or- 
thogonality of the eigenvectors A^ gives the stroboscopic 
propagator XJin{NTz) for N periods as 



UmiNTz 



Q 



-iEqNTz 



(11) 



For any integer N the system state is given by 
l^iNTz)) = U,„(iVTz)|*(0)). The similarity of 

Eq. (fTO| for a single-period propagator and Eq. ([TT|) for 
any A^ > 1 is a highly nontrivial fact demonstrating that 
Uin{NTz) — Uj^iTz) can be simply expressed by the 
right-hand-side in Eg.pTj). The stroboscopic approach 
allows us to study very accurately the long-time evolu- 
tion since the A^-period propagator pT|) is constructed 
explicitly in a finite algebraic form. Although this prop- 
agator describes the dynamics exactly, it allows to watch 
only the stroboscopic evolution rather than the entire 
continuous one. However, if we arc interested in slowly 
evolving phenomena such as chaos development [35| and 
Rabi oscillations which occur here on many periods of 
the driving field, the stroboscopic approach is fully jus- 
tified and highly efficient. The experiment [2| uses stro- 
boscopic approach with the intervals on the order of 100 
ns to measure the slow dynamics of the driven electron 
spin. 

The results of calculations of electron displacement at 
discrete times NTz are presented in Fig|3l As one can 
see in Fig 131 the time dependence of the displacement be- 
comes strongly nonperiodic with the typical values being 
considerably less than d. That is in a strong electric field 
the electron probability density redistribution between 
the dots is not complete. The motion can be qualita- 
tively analyzed in the pscudospin model of the charge 
dynamics, [2^ where the tunneling splitting is described 
by the Uz matrix, and the driving field is coupled to the 
CTa;— matrix. The decrease in the electron displacement 
with the increase in the electric field can be viewed as 
a suppressed spin precession in a strong periodic field or 
as a coherent destruction of tunneling. |32| This nonperi- 
odic behavior and decrease in the displacement eventu- 
ally result in a less efficient spin driving. It should be 
mentioned that the fast oscillations in Fig l2] which are in 
general absent in Fig 121 reflect the difference between the 
continuous time scale in the former Figure and the stro- 
boscopic Floquet times NTz in the latter one. Figure 
l3l clearly illustrates the role of the spin in the orbital dy- 
namics: the curves in FiglJJa) andlSlJb) are very different. 
Tracking of the system at stroboscopic times NTz raay 
not allow seeing the complete fast orbital dynamics, thus 




FIG. 3: (Color online) Stroboscopic time dependence of {x)/d 
for field / = 0.35 in a given time window: (a) B = 1.73 T, 
(b) B — 6.93 T. Solid lines serve only as a guide for the eye. 
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FIG. 4: (Color online) Stroboscopic time dependence of 
{(Jz{NTz)) for different external driving fields / marked near 
the plots for two different magnetic fields: (a) B = 1.73 T, 
vertical dashed line marks the operational definition of the 
spin-flip period, (b) same as in (a) for B = 6.93 T. 



masking some details. As a result, there is no simple way 
to describe this stroboscopic picture directly in terms of 
the Hamiltonian parameters. 

The slow long-term spin dynamics is presented in 
FigHl Here the "unit of time" Tz is short enough and 
the time dependence of (ctz) is accurately described by 
the stroboscopic approach. Since the spin dynamics is 
not strictly periodical and full spin flips do not always 
appear in this system, we use the operational definition 
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FIG. 5: (Color online) Nonlinear dependence of the Rabi spin- 
flip frequency on the electric field amplitude in multilevel sys- 
tem for two different magnetic fields (a) B = f .73 T, the char- 
acteristic spin-flip rate of the order of 50 MHz, (b) B = 6.93 
T, the characteristic spin-flip rate of the order of 400 MHz. 



of the "spin- flip" time Tg^: spin flip oecurs when spin 
component shows a broad minimum albeit accompanied 
by fast oscillations (see in FigUJ^a)). The fast dynamics 
in the spin- flip doublet shown in Fig |4] becomes slow with 
the field increase as a result of a weaker effective coupling 
of the states with different parity. The resulting spin be- 
havior, arising solely due to the SOC, is shown in FigU) 
The Rabi frequency for the spin-flip is smaller for some 
higher values of / (which we vary through Fig|31-Fig|3]) 
than for some weaker values of / in contrast to what 
can be expected for the weak fields employed, e.g. in 
the experiments [2| , being a manifestation of the gener- 
ally nonmonotonous behavior of the Rabi frequency on 
the electric field amplitude. In addition, in contrast to 
the simple Rabi oscillations, the flips become incomplete, 
with {(7z(t)) = — 1 never reached. These two qualitative 
effects are the results of the enhanced electron tunnel- 
ing between the potential minima: the spin precession in 
the driven interminima motion establishes corresponding 
spin dynamics and prevents the electric field to flip the 
spin efficiently. This effect makes a qualitative differ- 
ence to the model of Ref.Q, where electron is assumed 
to be always located in the orbital ground state near the 
minimum of the potential formed by the parabolic con- 
finement and weak external electric field. 

To present a broader outlook onto the dependence of 
the spin flip rate on the driving field, we plot in Fig|5]the 
spin fiip rate for B = 1.73 T and B = 6.93 T. In con- 
trast to the linear dependence for a conventional two-level 
Rabi resonance formula, one can see a strongly differ- 
ent much more complicated non-monotonous dependence 
in a multi-level structure, especially at high fields. The 
regime in FiglSJa) shows more irregularities since all four 
lowest levels are equidistant (Figl^l^a)) and involved in 



the resonance while in FiglSlJb) more regular dependence 
is observed, refiecting a simpler nature of the resonances 
here. 

We would like to mention here that the observed slow- 
ing down of spin dynamics can be seen on a more gen- 
eral ground, not restricted to the exact form of Eq.(|8l), 
as the Zeno effect of freezing evolution of a measured 
quantity. [3g-l33 Indeed, the operator —iaid/dx makes 
the orbital dynamics spin-dependent, and, as a result, 
performs the measurement of the Ci component [39|, |40| 
in the sense of von Neumann procedure. This can be seen 
in the evolution of a two-component wave function: [39| 

e-"*'^^^/^^</'(x)(Ci|l)+C-i|-l)) = 

(bix - at)Ci |1) + 0(a; + af)C-i |-1) , (12) 

where we took i — z as an example, Ci a-nd C-i corre- 
spond to ±1 eigenvalues of cr^, respectively, and a is the 
coupling constant. The SOC thus The SOC thus entan- 
gles the orbital and spin motion, destroys the coherent 
superposition of spin-up and spin-down states, and per- 
forms the von Neumann-like spin measurement by map- 
ping spin state on the electron position. This von Neu- 
mann measurement, is, however, different from the exper- 
imental measurement procedure applied, e.g. in Ref.[2|. 
The spin-orbit coupling coupling drives the coherent su- 
perposition of different spin components and at the same 
time, by constant strong measurement, destroys it lead- 
ing to a slow spin dynamics. 



CONCLUSIONS 



We have considered the interplay between the tun- 
neling and spin-orbit coupling in a driven by an exter- 
nal electric field one-dimensional single-electron double 
quantum dot. In the regime of the electric dipole spin res- 
onance, where the electric field frequency exactly matches 
the Zeeman transition, the complex interplay of these 
mechanisms results in two unexpected effects. The first 
effect is the nonmonotonous change in the Rabi spin os- 
cillations frequency with the electric field amplitude. The 
Rabi oscillations become much slower than expected for 
a two-level system. The second effect is the incomplete 
Rabi spin flips. This behavior results from the fact that 
the interminima motion establishes a competing spin dy- 
namics, leading to the physics somewhat similar to the 
Zeno effect, preventing a fast change in a measured quan- 
tity. These results indicating the slowdown and nonlin- 
earity of the spin resonance in multilevel systems can be 
useful for pointing out certain fundamental challenges for 
the future experimental and spintronics device applica- 
tions of phenomena based on spins in double quantum 
dots. 
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